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Abstract 

Here the effect of rotation up to third order in the angular velocity of a star on the p, f 
and g modes is investigated. To do this, the third-order perturbation formalism presented 
by Soufi et al. (1998) and revised by Karami (2008), was used. I quantify by numerical 
calculations the effect of rotation on the oscillation frequencies of a uniformly rotating (3- 
Cephei star with 12 Mq. For an equatorial velocity of 90 kms~^, it is found that the second- 
and third-order corrections for {l,m) = (5,-4), for instance, are of order of 0.07% of the 
frequency for radial order n = —3 and reaches up to 0.6% for n = —20. 
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1 introduction 

Pulsating stars on the upper main sequence, and particularly 6 Scuti and /3 Cephei stars, are 
rapid rotators as well as being multimode pulsators. The ratio e = Q/uj of the rotation rate 
to the typical frequency of oscillations co seen in these stars is no longer a small quantity as it is 
for e.g. the Sun. These stars have typically equatorial velocities ~ 100 kms~^, and oscillation 
periods from half to a few hours which implies e ~ 0.1, whereas for the Sun e ~ 10"'^. In order 
to achieve the full potential of asteroseismology for testing of models for upper main sequence 
stars a more careful treatment of the effect of rotation on oscillation frequencies is required. 

Rotation not only modifies the structure of the star but also changes the frequencies of 
normal modes. It removes mode degeneracy creating multiplets of modes. If the rotational 
angular velocity, Q, does not have any latitudinal dependence, and the rotation is sufficiently 
slow, the multiplets show a Zeeman-like equidistant structure. At faster rotation rates non- 
negligible quadratic effects in cause the position of the centroid frequency of multiplets to 
shift with respect to that of a non-rotating model of the same star. See Karami et al. (2003). 

Dziembowski & Goode (1992) derived a formalism for calculating the effect of differential 
rotation on normal modes of rotating stars up to second order. Soufi et al. (1998) extended 
the formalism of Dziembowski & Goode (1992) up to third order for a rotation profile that is a 
function of radius only. Soufi et al. (1998) found that near-degenerate coupling due to rotation 
only occurs between modes with either the same degree I (and different radial orders) or with 
modes which differ in degree by 2. In general they showed that the total coupling comes from 
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three distinct contributions: the Coriohs contribution, the non-spherically-symmetric distortion, 
and a coupling term which involves a combination of these two effects. 

Result of calculation of frequency corrections up to third order were presented for models 
of (5-Scuti stars by Goupil et al. (2001), Goupil & Talon (2002), Pamyatnykh (2003), and 
Goupil et al. (2004). Daszyhska-Daszkiewicz et al. (2002) studied the effects of mode coupling 
due to rotation on photometric parameters (amplitude and phase) of stellar pulsations. They 
reconfirmed the conclusion of Soufi et al. (1998) that the most important effect of rotation is 
coupling between close frequency modes of spherical harmonic degree, I, differing by 2 and of 
the same azimuthal order, m. 

Reese et al. (2006) studied the effects of rotation due to both the Coriolis and centrifugal 
accelerations on pulsations of rapidly rotating stars by a non-perturbative method. They showed 
that the main differences between complete and perturbative calculations come essentially from 
the centrifugal distortion. Suarez et al. (2006) obtained the oscillation frequencies include 
corrections for rotation up to second order in the rotation rate for 6 Scuti star models. Karami 
(2008) revised the third-order perturbation formalism presented by Soufi et al. (1998) because of 
some misprints and missing terms in some of their equations. Karami (2008) by the help of the 
revised formalism, calculated the effect of rotation up to third order on the oscillation frequencies 
of a uniformly rotating zero-age main-sequence star with 12 Mq. He concluded that for an 
equatorial velocity of 100 kms"^, the second- and third-order corrections for {l,m) = (2,2), for 
instance, are of order of 0.01% of the frequency for radial order n = 6 and reaches up to 0.5% 
for n = 14. 

In this paper, I use the third-order perturbation formalism according to Soufi et al. (1998) 
and revised by Karami (2008), hereafter Paper I. I carry out numerical calculations for the 
frequency corrections for a /? Cephei star with mass M = 12 Mq, Mq being the solar mass. The 
numerical results are presented in section 2. Section 3 is devoted to concluding remarks. 

2 Oscillations of a rapidly rotating P Cephei star 

In order to calculate the effect of rotation on normal modes, I consider a uniformly rotating, 
12 Mq, P Cephei model generated by the evolution code of Christensen-Dalsgaard (1982) (see 
also Christensen-Dalsgaard Sz Thompson 1999). The parameters of the model are listed at Table 
1. The value of central hydrogen abundance and existence of small convective core show that the 
model should be a quite evolved /3-Cephei star. The behavior of some of equilibrium quantities 
of the model against fractional radius, x = r/R are represented in Fig. 1; It shows that: 1) Close 
to the center up to radius x = 0.1, the star is in a convective regime where squared buoyancy 
frequency A'^^ < and, outside of this radius is in a radiative regime where iV^ > 0. There 
is a sharp peak in N'^ at x ~ 0.15. This happens because on one hand N'^ oc and on the 
other hand since the model has a very small convective core, .Rconv = 0.1 i?, hence the gravity, 
ge, increases sharply. 2) Spherically symmetric density p decreases smoothly from its maximum 
value to nearly zero near the surface at x ~ 0.4; 3) The absolute value of the non-spherically- 
symmetric correction to the density P22, see Eq. (10) in Paper I, does show a sharp peak at 
X ~ 0.15; 4) The absolute value of the non-spherically-symmetric correction to the gravitational 
potential (f>22, see Eq. (12) in Paper I, increases smoothly to its maximum value at the surface. 

2.1 Eigenfunctions 

The zero-order eigenfunctions arc computed from the zero order eigenvalue problem with the 
pulsation code of Christensen-Dalsgaard (see Christensen-Dalsgaard & Berthomieu 1991), mod- 
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ified according to Eqs. (20) to (24) in Paper I. 

In Fig. 2, the radial (y) and horizontal (z) components of the zero-order poloidal eigen- 
functions as well as rp^^'^S,r/{R'^Pc^'^) related to the radial energy density, where = ry is the 
radial displacement, are plotted against the fractional radius x = r/R for the selected modes 
with (/,m)=(5,-4) and n=(-7,...,-ll). The modes with n = (—8,-10) are the pure g-modes 
and with n = (—7, —9, —11) are the mixed g-modes. For the pure g-modes, the oscillations are 
mostly trapped near the center and also the horizontal amplitude of oscillations are comparable 
against the radial component. The horizontal component of buoyancy force which generates the 
horizontal amplitude, plays a important rule during the oscillation of a pure g-mode. At mathe- 
matical point of view a pure g-mode is mostly derived from a vector potential that its horizontal 
component has essential contribution in contrast with the corresponding radial component. This 
feature is clear particularly for polytropic models. See Sobouti (1980) and Sobouti & Rezania 
(2001). In the case of mixed g-modes, the oscillations are trapped between near the center and 
the middle part of the star. However the amplitude near the surface is decayed exponentially. 
Figure 2 also shows that: 1) The rapid oscillations in the pure and the mixed g-modes, occur at 
X = 0.1 due to existence of sharp peak in squared buoyancy frequency A^^ (see also Fig. 1). 2) 
Close to the center up to radius x = 0.1, where the star is in convective regime (iV^ < 0), the 
amplitudes of modes are decayed exponentially (see again Fig. 1). 

2.2 Eigenfrequencies and corrections 

The zero-order eigenfrequency, cto, is derived from numerical solutions of Eqs. (20) to (24) in 
Paper I by the modified pulsation code; note that using the eigensystem in Eqs. (20) to (24) the 
first-order frequency correction, cji, is implicitly included in uq (see Eqs. (15) and (16) in Paper 
I). The second- and third-order Coriolis contributions, (aj, crj), the second- and third-order 
non-spherically-symmetric distortions, (uP, 0-3*), and the third-order distortion and Coriolis 
coupling, a^, are derived from numerical integrations of Eqs. (52), (53), (64), (65), and (69) in 
Paper I. 

In Tables 2 to 9 the results of different contributions of frequency corrections due to effect of 
rotation up to third order are tabulated. In each table the selected p, f and g modes with (Z, m) = 
(5, -4), (2, 0), (2, 1), (2, 2), (2, -1), (5, 0), (5, 1), and (5, 4) with n = (-3, . . . , -20), (1, . . . , -5), 
(1, . . . , -5), (1, . . . , -4), (1, . . . , -5), (-2, . . . , -11), (-1, . . . , -10), and (0, . . . , -7), respectively, 
are considered. The modes with n > 1, n = 0, and n < are labelled by (pi, . . . ,Pn), f, and 
(5i,->ffn), respectively. 

Tables 2 to 9 show that: 1) The values of zero order eigenfrequency, ctq, and total frequency, 
fjtot, decrease when the radial mode number, n, decreases. 2) The order of magnitudes of {^2 ^ 
a^) are smaller than {aj , <tJ) by a factor of 10~^ to 10~^. Therefore one can concludes that the 
effect of Coriolis forces is dominant with respect to the centrifugal forces. 3) With increasing n, 
the frequency correction due to the distortion and the Coriolis coupling increases and decreases 
alternatively. 4) For the case of m = 0, there is no first and third-order frequency corrections. 

In Table 10, the results of third order frequency corrections for the case of two near degenerate 
modes, derived from Eqs. (75)-(76), are tabulated. The coupling exists only for the two near 
degenerate modes belonging to the same m but I differing by ±2. However there is no any 
selection rule for n. 

Note that in the numerical calculations, there are two substantial differences between the 
equations used in Soufi et al. (1998) and Paper I. In the formulation of Soufi et al. (1998), 

the density derivatives are eliminated through an integration by parts and the resulting surface 
terms are ignored. The surface terms become significant for higher-order modes, particularly in 
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the present model whose atmosphere is relatively thin. The other important difference which 
should be noted is that in Soufi et al. (1998) for computing the third-order correction terms the 
approximation z ~ Ut/CaQ, which is valid for the non-rotating case, is used. In Paper I, on the 
other hand, the exact relation Eq. (24) is used. The approximation of neglecting the second 
term in Eq. (24) everywhere, particularly near the surface, is not valid. The magnitude of this 
difference between the two approaches is more significant than the magnitude of the difference 
due to of the surface terms for the present model. If we include the surface terms in Soufi et 
al. (1998) and use the approximation z ~ yt/Cag in Paper I, the results of the two numerical 
approaches are in good agreement. 

3 Concluding remarks 

The third-order effect of rotation on the p, f and g modes for a uniformly rotating /3-Cephei 
star of mass 12 Mq has been investigated. The third-order perturbation formalism presented by 
Soufi et al. (1998) and revised by Karami (2008) was used. The zero-order eigenvalue problem 
was solved by pulsation code modified in this manner. Numerical calculations of oscillation 
frequencies were carried out for our selected model and second- and third-order frequency cor- 
rections due to Coriolis, non-spherically-symmetric distortion and Coriolis-distortion coupling 
were computed. For the case of m = 0, there is no first and third-order frequency corrections. 
Coupling only occurs between two poloidal modes with the same m and with / differing by or 2. 
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Table 1: Stellar parameters of a rotating /3-Cepliei star in solar units. M, Mconv, R, Rconv, Pc, 
and Pc are the total mass, the mass of convective core, the radius, the radius of convective core, 
the central pressure and density, and © denotes solar values; and e are the dimensionless 
mean angular velocity and the perturbational expansion coefficient; Tdyn, Ti-ot, and Vrot are 
the dynamical time scale (free fall time), the equatorial period and velocity, respectively. For 
comparison note that TQdyn = 0-5 h, Tqioi = 25 d, V0rot = 2 kms^^. is central helium 
abundance^ 



M = 12 Mq 
R = 8.92 Rq 
Pc = 2.85 X 10"^ Pc0 

= 1.7 8 X IQ-^ 
Tdyn = ^/RF/GM = 3.41 h 
y^^t = m = Ra^/Tdyn = 90 km/sec 



e = 



Pc = 5.39 X 10-2 pc0 
n/u = a^jl-K = 2.84 X 10-2 
Trot = 27ri?/Kot = 4.94 d 



Mconv = 0.18 M 




Xc = 0.3 
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P-Cephei Star (M=12 Solar Mass) 
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Figure 1: Dimensionless equilibrium quantities including squared buoyancy frequency A^^ = 
{gc/r)A with A = (l/ri)(dlnp/(ilnx) — (dlnp/dlnx), spherically and non-spherically- 
symmetric contributions p and P22 to the density, and non-spherically-symmetric gravitational- 
potential contribution (f)22, in units of GM/B?, pc and R'^O,'^, respectively, against fractional 
radius x = r/R for a /3-Cephei star model with M = 12Mq. 
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P Cephei Star a,m)=(5,-4) 
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Figure 2: Zero-order radial y (left) and horizontal z components (middle) as well as 
rp^^^Cr/iR'^pl^^), related to the energy density (right), against fractional radius x = r/R for 
selected g-modes with (Z,m) = (5, —4) and n = (—7, —8, —9, —10, —11) for the model described 
in Table 1. 
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Table 2: Values of eigenfrequency cq (including the 0{^}) contribution from rotation), second- 
and third-order Coriolis contributions aj and aj, second- and third-order non-spherically- 
symmetric distortions 02 and a^, third-order distortion and Coriolis coupling a^, and total 
frequency citot = Co + c^c corrected up to third order, for g modes in a /3-Cephei star with 
M = 12Mq, (/, m) = (5, -4) and n = (-3, -20). Here = crj + aj + + af + a§ is the 
total third-order frequency correction. All frequencies are in units of \/GM/ K'^ = 8.16 x lQ~^Hz. 
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Table 3: Same as Table 2, for p, f and g modes "with (i, m) = (2, 0) and n = (1, . . . , —5). There 
is no crj, (T3*, and (73" = for m = (see Eqs. (53), (65), and (69) in Paper I). 
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Table 4: Same as Table 2, for p, f and g modes with {I, m) = (2, 1) and n = (1, . . . , —5). 
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Table 5: Same as Table 2, for p, f and g modes with (/, m) = (2, 2) and n = (1, . . . , —4). 
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Pi 


1 


4.7531 


1.5883x10- 


~3 


-6.2091x10" 


-4 


-1.1783x10" 


-4 


3.8294x10" 


-5 


3.9523x10" 


~4 


4.7544 


f 





4.0472 


2.1026x10" 


-3 


-1.2523x10" 


-3 


-1.8862x10" 


-4 


9.3932x10" 


-5 


3.2290x10" 


-4 


4.0483 


91 


-1 


3.7510 


2.2968x10" 


-3 


-4.5920x10" 


-4 


-2.2129x10" 


-4 


3.6960x10" 


-5 


7.1757x10" 


-5 


3.7527 


92 


-2 


3.0857 


2.6160x10" 


-3 


-1.9895x10" 


-3 


-2.9955x10" 


-4 


1.8946x10" 


-4 


2.3278x10" 


-4 


3.0864 


93 


-3 


2.5394 


3.4410x10" 


-3 


-1.2564x10" 


-4 


-4.8923x10" 


-4 


1.4920x10" 


-5 


6.5731x10" 


-6 


2.5422 


9i 


-4 


1.9447 


5.0666x10" 


-3 


-2.3100x10 


-4 


-9.7698x10" 


-4 


3.7477x10" 


-5 


8.6389x10" 


-6 


1.9486 



Table 6: 


Same as 


Table 2, for 


P, 


f and g modes with {l,m 


) = 


(2, —1) and n 


= (1,...,- 


-5). 




Mode 


n 


CO 








^3^ 






-3^ 




Ctot 


Pi 


1 


5.2102 


4.4006x10" 


3 


3.6153x10-^ 


2.0474x10" 


-4 


1.1001x10-4 


3.5096x10" 


-5 


5.2186 


f 





4.5176 


5.1223x10" 


3 


2.3647x10"^ 


2.7408x10" 


-4 


8.4565x10--^ 


2.8493x10" 


-5 


4.5255 


9i 


-1 


4.2071 


5.2688x10" 


3 


7.9315x10-4 


3.0026x10" 


-4 


2.9067x10""^ 


1.7027x10" 


-5 


4.2135 


92 


-2 


3.5517 


6.5049x10" 


3 


1.3879x10-3 


4.4055x10" 


-4 


6.3833x10-5 


3.9379x10" 


-5 


3.5601 


93 


-3 


2.9898 


7.1166x10" 


3 


8.6542x10-'^ 


5.6663x10" 


-4 


4.3428x10-6 


1.0038x10- 


-5 


2.9976 


.94 


-4 


2.4105 


8.8416x10- 


3 


1.6696x10-'' 


8.7025x10- 


-4 


1.0610x10-^ 


4.4664x10- 


-5 


2.4204 


9r, 


-5 


2.104.") 


9.8081x10" 




1.9592x10-""' 


1.1018x10- 




1.3914x10-'* 


(i.()382xl0- 


-(i 


2.11.54 



Table 7: Same as Table 3, for g modes with {I, m) = (5, 0) and n = (—2, . . . , —11). 



Mode 


n 




-2^ 








Ctot 


92 


-2 


4.6430 


6.5692x10" 


-3 


5.1648x10" 


-3 


4.6548 


93 


-3 


4.1745 


7.4434x10" 


-3 


4.6698x10" 


-5 


4.1820 


94 


-4 


3.7054 


8.6005x10" 


-3 


1.0041x10" 


-3 


3.7150 


95 


-5 


3.4917 


8.5212x10" 


-3 


2.4579x10" 


-3 


3.5027 


96 


-6 


3.0567 


1.0158x10" 


-2 


3.1727x10" 


-5 


3.0669 


97 


-7 


2.8068 


1.1320x10" 


-2 


2.2605x10" 


-4 


2.8184 


98 


-8 


2.3747 


1.3073x10" 


-2 


2.4209x10" 


-5 


2.3878 


99 


-9 


2.2732 


1.3929x10" 


-2 


1.6712x10" 


-4 


2.2873 


910 


-10 


1.9343 


1.6049x10" 


-2 


1.9849x10" 


-5 


1.9504 


911 


-11 


1.8916 


1.6674x10" 


-2 


1.3491x10" 


-4 


1.9085 
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Table 8: Same as Table 2, for g modes with {I, m) = (5, 1) and n = (—1, . . . , —10). 



Mode 


n 








<^2^ 




-3^ 




-3^ 




-3^ 




ftot 


91 


-1 


5.2735 


5.7438x10- 


-3 


7.7440x10" 


-5 


-2.0385x10" 


-4 


-2.5593x10- 


-6 


-7.0613x10" 


-7 


5.2791 


92 


-2 


4.4767 


6.5970x10- 


-3 


2.6340x10" 


~3 


-2.6521x10" 


-4 


-9.7601x10- 


-5 


-1.3530x10" 


-4 


4.4855 


93 


-3 


4.0016 


7.5499x10- 


-3 


3.9642x10" 


-5 


-3.5099x10" 


-4 


-1.7135x10- 


-6 


-8.7762x10" 


-7 


4.0088 


54 


-4 


3.5245 


8.7660x10- 


-3 


8.3310x10" 


-4 


-4.7788x10" 


-4 


-4.2566x10- 


-5 


-3.4359x10" 


-5 


3.5335 


.95 


-5 


3.3328 


8.6929x10- 


-3 


1.6871x10" 


-3 


-4.5513x10" 


-4 


-8.0691x10- 


-5 


-1.0396x10" 


-4 


3.3426 


96 


-6 


2.8839 


1.0492x10- 


-2 


2.6152x10" 


-5 


-6.7611x10" 


-4 


-1.5664x10- 


-6 


-9.6720x10" 


-7 


2.8938 


97 


-7 


2.6261 


1.1821x10- 


-2 


1.5839x10" 


-4 


-8.6786x10" 


-4 


-1.0907x10- 


-5 


-7.8477x10" 


-6 


2.6372 


98 


-8 


2.2021 


1.3773x10- 


-2 


1.9341x10" 


-5 


-1.1620x10" 


-3 


-1.5166x10- 


-6 


-9.7065x10" 


-7 


2.2147 


99 


-9 


2.0938 


1.4799x10- 


-2 


1.0486x10" 


-4 


-1.3547x10" 


-3 


-8.9901x10- 


-6 


-7.3938x10" 


-6 


2.1073 


910 


-10 


1.7617 


1.7257x10- 


-2 


1.5214x10" 


-5 


-1.8197x10" 


-3 


-1.4909x10- 


-6 


-9.8237x10" 


-7 


1.7771 
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Table 9: Same as Table 2, for f and g modes with {I, m) = (5, 4) and n = (0, . . . , —7). 



Mode 


n 












-3^ 




-3^ 




-3^ 




t^tot 


f 





4.9968 


2.6949x10- 


-3 


2.9857x10" 


-3 


-3.8428x10" 


-4 


-3.9355x10- 


-4 


2.5204x10- 


-4 


5.0020 


91 


-1 


4.7510 


3.0994x10" 


-3 


-4.5756x10" 


-5 


-4.8592x10" 


-4 


6.7139x10- 


-6 


-7.5470x10" 


-7 


4.7535 


92 


-2 


3.9831 


3.3532x10 


-3 


2.5387x10" 


-3 


-5.9624x10 


-4 


-4.1671x10- 


-4 


3.0860x10- 


-4 


3.9883 


93 


-3 


3.4828 


4.1984x10" 


-3 


-2.1823x10" 


-5 


-8.9241x10" 


-4 


4.3348x10- 


-6 


-9.4382x10" 


-7 


3.4861 


94 


-4 


2.9964 


4.8000x10" 


-3 


-1.1922x10" 


-4 


-1.1736x10" 


-3 


2.7171x10- 


-5 


1.2330x10- 


-4 


3.0000 


95 


-5 


2.8464 


4.8577x10" 


-3 


-1.6612x10" 


-4 


-1.2257x10" 


-3 


3.8867x10- 


-5 


1.0211x10- 


-4 


2.8500 


96 


-6 


2.3658 


6.2203x10" 


-3 


-1.2587x10" 


-5 


-1.9446x10" 


-3 


3.6761x10- 


-6 


-9.7763x10" 


-7 


2.3701 


97 


-7 


2.0828 


7.6637x10" 


-3 


-1.8813x10" 


-5 


-2.8324x10" 


-3 


6.5530x10- 


-6 


-9.7497x10" 


-6 


2.0876 



Table 10: Values of zero-order cigcnfrequency do, total frequency a± corrected up to third 
order, total third-order frequency correction Aa± = {a± — o"o) due to rotation and coupling, 
expansion coefficients ^^^^ and A'^^ normalized to 1, i.e., Al + A2 = 1, for selected pairs of 
near-degenerate poloidal modes in a /3-Cephei star with M = 12Mq. All frequencies are in units 
of s/GMfW = 8.16 X lO^^Hz. 



m. 


/ 


n 


coupling 






Acr± 


















3 


P3 


5.1186 


5.1271 


8.4915x10" 


-3 


9.98528x10- 


-1 


5.42353x10- 


-2 





2 


1 


Pi 


5.0537 


5.0640 


1.0282x10" 


-2 


5.72211x10- 


-2 


-9.98362x10" 


-1 





2 


-2 


92 


3.3908 


3.4016 


1.0888x10" 


-■2 


9.99999x10- 


-1 


1.65116x10- 


-3 





4 


-4 


9i 


3.4704 


3.4795 


9.1303x10" 


-3 


8.92439x10- 


-3 


-9.99960x10" 


-1 


1 


3 


-3 


93 


3.3049 


3.3141 


9.2414x10" 


-3 


9.98573x10- 


^1 


5.34025x10- 


-2 


1 


5 


-5 


95 


3.3328 


3.3426 


9.8344x10" 


-3 


4.56655x10- 


-2 


-9.98957x10" 


-1 


2 


3 


-3 


93 


3.1457 


3.1511 


5.3915x10" 


-3 


9.98319x10- 


-1 


5.79619x10- 


-2 


2 


5 


-5 


95 


3.1729 


3.1809 


8.0497x10" 


-3 


5.62952x10- 


-2 


-9.98414x10" 


-1 


3 


3 


-3 


93 


2.9905 


2.9908 


3.3252x10" 


-4 


9.98752x10- 


-1 


4.99354x10- 


-2 


3 


5 


-5 


95 


3.0111 


3.0171 


5.9922x10" 


-3 


5.81273x10- 


-2 


-9.98309x10" 


-1 
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